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t-H ■ ABSTRACT 

> ■ 

Aims. The aim of the present work is the construction of a mass-selected galaxy cluster sample based on weak gravitational 
lensing methods. This sample will be subject to spectroscopic follow-up observations. 

Methods. We apply the mass aperture statistics and a derivative of it to 19 square degrees of high quality, single colour wide 
field imaging data obtained with the WFIOMPG/ESO 2.2m telescope. For the statistics a family of filter functions is used that 
approximates the expected tangential radial shear profile and thus allows for the efficient detection of mass concentrations. 
Results. We identify 158 possible mass concentrations. This is the first time that such a large and blindly selected sample is 
published. 72 of the detections are associated with concentrations of bright galaxies. For about 22 of those we found spectra in 
the literature, indicating or proving that the galaxies seen are indeed spatially concentrated. 15 of those were previously known 
to be clusters or have meanwhile been secured as such. We currently follow-up a larger number of them spectroscopically to 
Q ■ obtain deeper insight into their physical properties. The remaining 55% of the possible mass concentrations found are not 
S—j associated with any optical light, or could not be classified unambiguously. We show that those "dark" detections are to a 
CZ) , significant degree due to noise, and appear preferentially in shallow data. 
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1. Introduction peaks were reported so far (see Erben et al. 2000; Umetsu 

& Futamase 2000; Miyazaki et al. 2002; Dahle et al. 2003 

After the recent determination of the fundamental cosmo- Wittman et aL 2 001, 2003; Schirmer et al. 2003, 2004 

logical parameters (Spergel et al. 2006) a profound under- Hetterscheidt et al. 2005; von der Linden et al. 2005 

standing of the dark and luminous matter distribution in wittman et aL 2 005, for some examples), a very small 

the Universe is one of the key problems in modern cosmol- number compared to the more than 10000 clusters known 

ogy. Galaxy clusters are at the centre of attention in this to date (see Lopes et aL 20 04, for 9900 cluster candidates 

context since they indicate the largest density peaks of the Qn northern sky within z — 1 5) 
cosmic matter distribution. Their masses can be predicted 

by theory, and gravitation dominates their evolution. Thus The selection of mass concentrations using the shear 
clusters are prime targets for the comparison of observa- caused by their weak gravitational lensing effect suffers 
tion against theory. For this purpose they should preferen- from a number of disadvantages. The most important one 
tially be selected by their mass instead of their luminosity is the high amount of noise contributed by the intrinsic el- 
in order to avoid a selection bias against underluminous lipticities of lensed galaxies. This largely blurs the view of 
members. Weak gravitational lensing methods select clus- the cosmic density distribution, letting peaks disappear, 
ter candidates based solely on their mass properties, but and fakes peaks where there actually is no overdensity. It 
this method has not yet been used systematically on a can only be beaten down to some degree by deep observa- 
very large scale (some 100 or 1000 square degrees) due to tions in good seeing. The other disadvantage is that any 
a lack of suitable data. Only a few dozen shear-selected mass along the line of sight will contribute to the signal, 
giving rise to false peaks. Such projection effects or cos- 
Send offprint requests to: M. Schirmer mic shear can only be eliminated or at least recognised if 
* Based on observations made with ESO Telescopes at the redshift information of either the lensed galaxies, or the 
La Silla Observatory matter distribution in the field is available. Cosmic shear 
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can act as a source of noise (see for example Maturi et al. 
2005), but it has recently been shown by Maturi et al. 
(2006) that this can partly be filtered out. 

In the present work we use the aperture mass statistics 
(Mn P ) (Schneider 1996, hereafter S96) and a derivative of 
it for the shear-selection of density peaks, based on 19 
square degrees of sky coverage. The purpose of the work 
is to establish a suitable filter function for M ap , and then 
apply it to an (inhomogeneous) set of data. The sample 
returned is currently the largest sample of shear-selected 
cluster candidates, yet is dwarfed by the total number of 
galaxy clusters known. 

The outline of this paper is as follows. In Sect. 2 we give 
an overview of the data used, concentrating on its qual- 
ity and showing the usefulness for this analysis. Section 3 
contains a discussion of our implemented version of M ap , 
particularity with regard to the chosen spatial filter func- 
tions. We also introduce a new statistics, deduced from 
M ap . In Sect. 4 we present and discuss our detections, 
and we conclude in Sect. 5. 

Throughout the rest of this paper we use common weak 
lensing notations, and refer the reader to Bartelmann & 
Schneider (2001) and Schneider et al. (2006) for more de- 
tails and technical coverage. 

2. Observations, data reduction, data quality 

2.1. The survey data 

The GaBoDS was conducted with the wide field imager 
(hereafter WFI@2.2) of the 2.2m MPG/ESO telescope in 
the i?-band. It is to about 80% a virtual survey, which 
means most of the data was not taken by us but collected 
from the ESO archive. There are five main data sources: 

— The ESO Distant Cluster Survey, consisting of 12 
pointings that have been selected because z = 
0.5... 0.8 cluster candidates have been identified in 
them using colour techniques. These exposures are 
rather shallow and we are thus insensitive to such dis- 
tant objects (Fig. 8). Therefore, these fields do not 
impose a selection effect and qualify for our survey. 

— Our ASTRO VIRTEL 1 project for data mining the 
ESO archive, which yielded 30 random pointings. As 
compared to the rest of the fields, many of them have 
unsuitable PSF properties, e.g. due to less careful fo- 
cusing, so that we selected the 19 best fields. Details 
of this approach and its implementation are given in 
Schirmer et al. (2003) and Micol et al. (2004). 

— The COMBO-17 data makes for the deepest survey 
part. Four of the five fields (Sll, A901, SGP and FDF) 
contain known clusters. 

— The EIS Deep Public Survey, contributing 9 empty 
fields. 

— Our own observations of 13 empty fields. 

The final sky coverage is 18.6 square degrees, spread 
over 58 fields which are suitable for our weak gravitational 

1 http:/ /www.stecf.org/astrovirtel/ 



lensing analysis. The overall sky distribution of the fields 
is shown in Fig. 1, and further characteristics are sum- 
marised in Table B.l. 

2.2. Data reduction 

For the reduction of this specific data set we developed a 
stand-alone, fully automatised pipeline (THELI), which 
we made freely available to the community. It is capable 
of reducing almost any kind of optical, near-IR and mid- 
IR imaging data. A detailed description of this software 
package is found in Erben et al. (2005), in which we inves- 
tigate its performance on optical data. One of the main 
assets of this package is a very accurate astrometric solu- 
tion that does not introduce any artificial PSF distortions 
into the coadded images. 

The only difference between the current version of 
THELI and the one we used for the reduction of the sur- 
vey data over the last years, is that for the image coad- 
dition ElSdrizzle was used, which is meanwhile replaced 
by Swarp. The latter method leads to a 4% smaller PSF 
in the final image in case of superb intrinsic image seeing 
as in our survey. The PSF anisotropy patterns themselves 
are indistinguishable between the two coaddition methods. 
Since the natural seeing variations (Fig. 1) in our images 
are much larger than those 4%, our analysis remains un- 
affected. 

2.3. Image seeing and PSF properties 

The image seeing and the PSF anisotropy are critical for 
weak lensing measurements. They dilute and distort the 
signal and determine how well the shape of the lenscd 
galaxies can be recovered. As can be seen from Fig. 1, 
80% of our coadded images have sub-arcsecond image see- 
ing, and 20% are around l'.'O, with an average of 0'.'86. 
Thus, we reach the sub-arcsecond seeing regime very well 
which is commonly regarded as mandatory for weak lens- 
ing purposes. 

PSF anisotropics are rather small and usually well- 
behaved with WFI@2.2, which we have demonstrated 
several times (Schirmer et al. 2003, 2004; Erben et al. 
2005, for example). With a well-focussed telescope, 1% 
of anisotropy in sub-arcsecond seeing conditions can be 
achieved, with a long-term statistical mode of around 2% 
(see Fig. 3). Discontinuities in the PSF are largely ab- 
sent across chip borders. Slightly defocused exposures ex- 
hibit anisotropies of 3 — 5%. We rejected individual ex- 
posures from the coaddition if one or more CCDs had an 
anisotropy of larger than 6% in either of the cllipticity 
components E\ or e 2 - These anisotropies arise from astig- 
matism, and they flip by 90 degrees if one passes from an 
intrafocal to an extrafocal exposure (see Schirmer et al. 
2003, for an example of this, and Fig. B.2 for a typical 
PSF anisotropy pattern in our data). 

Since such intra- and extrafocal exposures are 
roughly equally numbered for a larger set of exposures, 
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Fig. 1. Left: Sky distribution of the GaBoDS fields. The size of the symbols indicates the covered sky area (not to scale), 
with one WFI shot covering 0.32 square degrees. All fields are at high galactic latitude. Right: Image seeing of the 58 coadded 
WFI@2.2 mosaics used for the lensing analysis. The average seeing is 0'.'86. 



anisotropies due to defocusing average out in the coadded 
images. Those have on average an anisotropy of 1 — 2% 
with a similar amount in the rms (variation of the PSF 
across the field). This is illustrated in Fig. 2, where we 
show the combined PSF anisotropy properties for all coad- 
ded mosaic images. The left panel shows the uncorrected 
mean ellipticity components, having (ei) = 0.0138 and 
(£2) = —0.0002. These anisotropies are small and become 
I (£1,2) I < 0.001 after PSF correction for all mosaics. 

The middle plot of Fig. 2 shows the uncorrected rms 
values for (£1.2), i.e. the deviations of the PSF from a con- 
stant anisotropy across the field. The rms of both compo- 
nents peaks around 1% and is reduced by a factor of 2 after 
PSF correction (right panel). The tail of the distribution 
seen in the middle plot essentially vanishes. 

To evaluate the remaining residuals from PSF cor- 
rection more quantitatively, we calculated the correlation 
function between stellar ellipticity and shear before and 
after PSF correction (Fig. 4), separately for the various 
survey data sources and over all galaxy positions, 

N 

^I'flia) = l/JV £ e *i P 2 Ol (<?h n ,2(0) . (1) 

i=l 

Here 71,2 are the PSF-corrected ellipticities of the galax- 
ies which serve as an unbiased estimator of the shear, and 
£*'2° are the components of the PSF correction polyno- 
mial at the position 8 of a galaxy. We find that the corre- 
lation between stellar PSF and shear is greatly reduced by 
the PSF correction, yet some non-zero residuals remain as 
we expected. Another test for residual systematics from 
PSF correction in the final lensing catalogues is the M ap 
(defined in equation (3)) two-point cross-correlation func- 
tion of uncorrected ellipticities of stars and corrected el- 
lipticities of galaxies (see Fig. 5), 

JVficlda 

« p Mf p al > = £ (M a * p M a f >, . (2) 



We find small non-zero residuals with different behaviour 
for the various survey parts. In particular, these residuals 
become increasingly non-zero with growing aperture size 
for our own observations and part of the ASTROVIRTEL 
fields. We will come back to these two points later during 
our analysis in Sect. 4, concluding that they do not affect 
our cluster detection method in a noticeable way. 

To summarise, our PSF correction effectively takes out 
mimicked coherent shear patterns from the data, yet small 
residuals remain, which are much smaller than the coher- 
ent shear signal (a few percent) we expect from a typical 
cluster at intermediate redshift range (z = 0.1 . . . 0.4). 

2.4. Catalogue creation: object detection and 
calculation of lensing parameters 

Our catalogue production can be split into three parts, the 
object detection, the calculation of the basic lensing quan- 
tities, and a final filtering of the catalogue obtained. An 
absolute calibration of the magnitudes with high accuracy 
is not needed for this work since we are mostly interested 
in the shapes of the lensed galaxies but not in their fluxes. 
We adopted photometric standard zeropoints which bring 
us, conservatively estimated, to within 0.05 — 0.1 mag of 
the real photometric zeropoint. 

For the detection process we use SExtractor (Bertin & 
Arnouts 1996) to create a primary source catalogue. The 
weight map created during the coaddition (see Erben et al. 
2005) is used in this first step, guaranteeing that the highly 
varying noise properties in the mosaic images are correctly 
taken into account. This leads to a very clean source cata- 
logue that is free from spurious detections. The number of 
connected pixels (DETECT_MINAREA) used in this work 
for the object detection was 5, and we set the detection 
threshold (DETECT.THRESH) to 2.5. These thresholds 
are rather generous. 10 — 25% of the objects detected are 
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Fig. 2. Left: The average £1,2 ellipticity components for the PSF of each mosaic, before PSF correction. The £2 component is 
symmetric around zero, whereas the £1 component scatters broader around a value of 1.4%. We do not show the corresponding 
plot after PSF correction, since all |ei,2| < 0.001. Middle: The scatter of the two components around their average values. This 
indicates the deviation from a constant PSF anisotropy across the field and the noise level. Right: Same as in the middle, but 
after applying our PSF correction scheme. As a result, not only does (£1,2) get largely removed, also their rms is significantly 
lowered. 



rejected again in the later filtering process since they were 
too small or too faint for a realiable shear measurement, 
or other problems appeared during the measurement of 
their shapes or positions. 

In the second step we calculate the basic lensing quan- 
tities using KSB (Kaiser et al. 1995). This includes PSF 
anisotropy correction, and the recalculation of the objects' 
first brightness moments since SExtractor yielded posi- 
tions with insufficient accuracy 2 for the purpose of our 
analysis. Details of our implementation of KSB are given 
in Erben et al. (2001), and a mathematical description of 
the PSF anisotropy correction process itself can be found 
in Bartclmann & Schneider (2001). 

2.5. Catalogue creation: filtering 

The third step consists of filtering the catalogue for various 
unwanted effects. To avoid objects in the catalogue that 
are in the immediate vicinity of bright stars, we explicitly 
set BACK.TYPE = MANUAL and BACK.VALUE = 0.0 
(our coaddcd images are sky-subtracted) in the SExtractor 
configuration file. Thus SExtractor is forced to assume 
a zero sky background and does not model the haloes 
around bright stars as sky background. This proved to 
be a very efficient way of automatically masking brighter 
stars and the regions immediately surrounding them (see 
e.g. the lower right panel of Fig. B.4). Further filtering 
on the SExtractor level is done by excluding all objects 
that are flagged with FLAG > 4 and those with negative 
half-light radii. 

On the KSB level we filter such that only galaxies for 
which no problems in the determination of centroids oc- 



2 This has been overcome in the recent SExtractor releases 
(v2.4.3 and higher), using a Gaussian weighting function in the 
process. 



curred remain in the catalogue. Galaxies with half-light 
radii (rjj smaller than 0.1 — 0.2 pixels than the left ridge 
of the stellar branch in an rh-mag-diagramme are rejected 
from the lensing catalogue, as are those with exceed- 
ingly bright magnitudes or a low detection significance 
(i/ max < 10). See the left panel of Fig. 6 for an illustra- 
tion of these cuts. From the same panel it can be seen 
that a significant number of galaxies have half-light radii 
comparable to or a bit smaller than the PSF, which makes 
their shape measurement noisier. Yet their number is large 
enough so that the shear selection of galaxy clusters profit 
significantly when these objects are included in the calcu- 
lation. By including these objects, we gain 10 — 25% in 
terms of the number density of galaxies, and 3 — 10% in 
terms of signal-to-noise of the detections. 

Furthermore, all galaxies with a PSF corrected modu- 
lus of the ellipticity larger than 1.5 are removed from the 
catalogue (the ellipticity can become larger than 1 due to 
the PSF correction factors, but is then downweighted) , as 
are those for which the correction factor (TrP s ) _1 > 5 
(see Erben et al. 2001). The fraction of rejected galaxies 
due to the cut-off in P g is relatively small, as can be seen 
from the right panel in Fig. 6. 

An overall impression of the remaining objects in 
the final catalogues is given in Fig. B.l. In total, typi- 
cally 10 — 25% of the objects are rejected from the ini- 
tial catalogue due to the KSB filtering steps. The re- 
maining average number density of galaxies per field is 
n <~ 11 arcmin~ 2 (min: 6, max: 28), not corrected for the 
S Extractor-masked areas (as described at the beginning of 
this section; on the order of 5% per field). For the width of 
the ellipticity distribution we measure o e — 0.34 for each 
of the two ellipticity components, averaged over all survey 
fields. Both n and o e determine the signal-to-noise of the 
various mass concentrations detected. 
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Fig. 3. Characteristic PSF anisotropics (|e*|) (modulus of the 
ellipticity of stars, averaged over the entire CCD) of 700 ran- 
domly chosen exposures for CCDs 1 and 8 of the WFIQ2.2 
detector mosaic. The remaining six CCDs have somewhat bet- 
ter properties. The exposure times were in the range of 300-500 
seconds, and the zenith distances were smaller than 40 degrees. 
The differences between the distributions arise from slightly 
different spatial orientations with respect to the focal plane 
(the CCDs are confined within ±20 microns), thus responding 
differently to the focus. Anisotropies larger than 1 — 2% are 
mostly due to an increasingly defocused instrument. The ex- 
posures used for this statistics were taken by a dozen different 
observers spread over more than a year, and give an idea of 
the quality of the archival data. 

3. Detection methods 

3.1. S-statistics and an optimal filter Q 

We base our detection method on the aperture mass statis- 
tics (M ap ) as introduced by S96. M ap can be written as a 
filtered integral of the tangential shear j t , 



M ap (0 o ) = f dV7t(^;0o) Q(<p). 



(3) 



Originally, the idea of M ap is to obtain a measure of the 
mass inside an aperture independent of the mass sheet 
degeneracy in the weak lensing case. Written in the form 
(3) we can also simply interpret it as the filtered amount 
of the tangential shear around a fiducial point 9q on the 
sky, where 7t(<£>; do) is the tangential shear at position <p 
relative to Oo, and Q is some radially symmetric spatial 
filter function. 

The variance of M ap for the unlensed case, respectively 
the weak lensing regime, is given by 
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Fig. 4. Shown is the correlation between stellar ellipticity and 
measured shear before (left panel) and after (right panel) PSF 
correction for the 58 survey fields. The median improvement is 
a factor of 3, but residuals are still present. 
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Fig. 5. M a p two-point cross-correlation function of uncor- 
rected ellipticities of stars and corrected ellipticities of galaxies, 
binned for the various survey data sources. The analysis of the 
ASTRO VIRTEL fields has been split up into the B8- and C0- 
fields, since these exhibit different properties. 



where a e is the dispersion of intrinsic source ellipticities, 
and n the number density of background galaxies. The 
integration is performed over a finite interval since we will 
have Q ~ for ■& > 9, i.e. for radii larger than the aperture 
size chosen. For the application to real data we will replace 
this integral in Sect. 3.2 by a sum over individual galaxies. 

We then define the S-statistics, or the S/N for M ap 
respectively the measured amount of tangential shear, as 



S(Mo) 



r o 9 d 2 7? 7t (fl;go)Q(fl) 



(5) 
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Fig. 6. Left panel: Stars appear as a vertical branch in a rh-mag plot. Those brighter than R — 16.5 saturate the detector and 
thus increase in size. The solid line encircles the galaxies which are used for the lensing analysis. The upper curved line indicates 
a cut in detection significance (f max > 10), which has proven to work significantly better than a constant cut on the faint end 
of the magnitudes. Right panel: c m „ against P s . Objects left of the indicated threshold are rejected from the lensing catalogue. 
Typically 1% of the galaxies are removed during this step. 



Here 8 is the aperture radius, and i? measures the distance 
inside this aperture from its centre at n . This expression 
gets a bit simplified due to assumed radial symmetry (see 
Appendix A). Hereafter, we will call the 2-dimensional 
graphical representation of the 5-statistics the S-map. If 
we plot the ^-statistics for a given mass concentration as 
a function of aperture size, then we refer to this curve as 
the S-profile. 

The filter function Q that maximizes S for a given 
density (or shear) profile of the lens can be derived us- 
ing either a variational principle (Schirmer 2004), or the 
Cauchy-Schwarz inequality (S96). It is obtained for 

Q(0)«7tW, (6) 

where 7t(i?) is the tangential shear of the lens, averaged 
over a circle of angular radius d. This is intuitively clear, 
since a signal (i.e., 7 t ) with a certain shape is best picked 
up by a similar filter function (Q). We present in Sect. 
3.4 a mathematically simple family of filter functions that 
effectively fulfills this criterion for the NFW mass profile. 

3.2. Formulation for discrete data fields 

For the application to real data, the previously introduced 
continuous formulation has to be discretised. First, we re- 
place the tangential shear j t with the tangential ellipticity 
e t , which is in the weak lensing case an unbiased estimator 



of 7 t . Thus we have 

1 N 

M a p = - £ " Qi , (7) 

2 — 1 

where n is the number density of galaxies, and N the total 
number of galaxies in the aperture. Introducing individual 
galaxy weights un as proposed by Erben et al. (2001), this 
becomes 

M ap = — V S U W t Qi , (8) 

where A is the aperture area, previously absorbed in the 
number density n. 

The noise of M ap can be estimated from M ap itself as 
was shown by Krusc & Schneider (1999) and S96. In the 
weak lensing case its variance evaluates as 

(7 2 (M ap ) = (Map 2 ) - (M ap ) 2 -(M a p 2 ), (9) 

since one expects (M ap ) = 0. Substituting with equation 
(7) yields 

<r 2 (M ap ) = \ J2(e u e tJ )Q l Q ] = \ J2(e tl 2 )Q l 2 

= 2^EN a «* 9 . (io) 

i 

where we used 

(et 2 ) = \ \e\ 2 , (11) 
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and the fact that e ti and e t j are mutually independent 
and thus average to zero for i =/= j. The Qi are constant 
for each galaxy and thus can be taken out of the averaging 
process. Again, in the case of individual galaxy weights 
this becomes 



! (M ap ) 



A 2 



2 (E< 



Q% 



(12) 



using the fact that the Wi are constant for each galaxy like 
the Qi. Therefore, we obtain for the ^-statistics 




S = 



V2 J2i £ t» Wj Qi 



(13) 



which grows like \/~N . 



3.3. Validity of the M ap concept for our data 

If M ap is evaluated close to the border of an image or on a 
data field with swiss-cheese topology due to the masking of 
bright stars, then the aperture covers an 'incomplete' data 

field. Therefore, the returned value of M ap does no longer 3,4. Filter functions 
give a result in the sense of its original definition in S96, 
which was a measure related to the filtered surface mass 
density inside the aperture. Yet it is still a valid measure 
of the tangential shear inside the aperture, including the 
5/Ar-estirriate in equation (13), and can thus be used for 
the detection of mass concentrations. 

Since the number density of background galaxies in- 
side an aperture is not a constant over the field due to the 
masking of brighter stars (and the presence of the field 
border), we have to check for possible unwanted effects. 
As long as the holes in the galaxy distribution are small 
compared to the aperture size, and as long as their number 
density is small enough so that no significant overlapping 
of holes takes place, the effects on the S-statistics are neg- 
ligible (see Fig. B.4). In fact, the decreased number density 
just leads to a lowered significance of the peaks detected 
in such areas, without introducing systematic effects. 

If the size of the holes becomes comparable to the aper- 
ture, spurious peaks appear in the 5*-map at the position 
of the holes. This is because the underlying galaxy pop- 
ulation changes significantly when the aperture is moved 
to a neighbouring grid point. When such affected areas 
were present in our data, then we excluded them from the 
statistics and masked them in the /S-maps, even though 
these spurious peaks are typically not very significant 
(<~ 2<t). Our threshold for not evaluating the 5-statistics 
at a given grid point is reached if the effective number 
density of the galaxies in the aperture affected is reduced 
by more than 50% due to the presence of holes (or the im- 
age border) . Spurious peaks become very noticeable if the 
holes cover about 80% of the aperture. This is rarely the 
case for our data unless the aperture size is rather small 
(2' ), or a particular star is very bright. We conclude that 
our final statistics is free from any such effects. 



Fig. 7. Left panel: Qnpw (solid line) , shown with some repre- 
sentations of Qtanh- The exponential cut-off E(x) for small 
and large radii is not introduced in this plot in order to 
show the differences between the two filter types better. As 
can be seen, the Qtanh filter is a very good approximation 
for Qnpw, giving slightly more weight to smaller radii. Right 
panel: Qtanh with the cut-off introduced by E(x) at both 
ends, plotted against the radial coordinate ■&. The dashed line 
(9 = 1) can be compared directly to the dash-dotted line in 
the left panel. 



As expressed in (6), an optimal filter function should re- 
semble the tangential shear profile. In the following such a 
filter Q(x) is constructed, assuming that the azimuthally 
averaged shear dependence is caused by an NFW density 
profile (Navarro et al. 1997) of the lensing mass concentra- 
tion. We set x := $/6, with x being the projected angular 
separation $ on the sky from the aperture centre, in units 
of the aperture radius 9. By varying 8, shear patterns re- 
spectively mass concentrations of different extent can be 
detected. 

Wright & Brainerd (2000) and Bartclmann (1996) de- 
rived an expression for the tangential shear of the univer- 
sal NFW profile. Based on their finding we can construct 
a new filter function Qnfw over the interval x G [0,1], 
having the shape 



Qnfw(x) 



4(3jy 2 -2) 



f -41n2 

4(3y 2 -2) 



y (y 



arctanhi 



arctan* 



y 2 1 1-y 2 



(14) 

(»<1) 
(» = 1) 

(y>i) 



Here we defined y :— x/x c , with x c being a dimcnsionless 
parameter changing the width (and thus the sharpness) 
of the filter over the interval x € [0, 1], in the sense that 
more weight is put to smaller radii for smaller values of 
x c 3 . This expression is smooth and continuous for y = 1, 
and approaches zero as ln(y)/y 2 for i/> 1. 

Due to the mathematical complexity of QnfW) the cal- 
culation of the S-statistics is rather time consuming for a 



Thus x c is in analogy to the NFW scale radius r s 
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field with <~ 10 4 galaxies. We introduce an approximat- 
ing filter function with simpler mathematical form that 
produces similarly good results as Qnfw- It is given by 

^ , ^ . . tanhfx/xc) 

Qtanh(x) = E(x) p-^, (15) 

X j x c 

having the 1 jx dependence of a singular isothermal sphere 
for x ^> x c . The hyperbolic function, having a oc x de- 
pendence for small x, absorbs the singularity at x = 
and approaches 1 for growing values of x. The pre- factor 
E{x) is a box filter, independent of x c , with exponentially 
smoothed edges. It reads 

E( X ) = l + e 6- 150x _|_ e -47 + 50x 

and lets Q drop to zero for the innermost and outermost 
10% of the aperture, while not affecting the large rest 
(see right panel of Fig. 7). It is introduced because both 
Qnfw and Qtanh (without this pre-factor) are not zero 
at the centre nor at the edge of the aperture. This cut-off 
is very similar to that introduced in [S96] for a differ- 
ent radial filter profile. It suppresses stronger fluctuations 
when galaxies enter (or leave) the aperture, receiving sig- 
nificant non-zero weight. It also avoids assigning a large 
weight to a few galaxies at the aperture centre, which as 
well can lead to significant fluctuations in the S'-statistics 
when evaluated on a grid. The effect of E(x) is rather mild 
though, since usually several hundred galaxies are covered 
by one aperture unless it is of very small size so that E(x) 
becomes important. 

We thus have a filter function based upon the two- 
dimensional parameter space (8 ,x c ). The differences be- 
tween Qtanh and Qnfw are indistinguishable in the noise 
once applied to real data, so that we do not consider Qnfw 
henceforth. 

It is not the first time that M ap filters following the 
tangential shear profile are proposed or used. We have 
already utilised the filter in equation (15) to confirm a se- 
ries of luminosity-selected galaxy clusters (Schirmer et al. 
2004). Before that, Padmanabhan et al. (2003) approxi- 
mated Qnfw with 

g PAD(x) = i^±£)__l__ 1/(1 + xf , (17) 

which was later-on modified by Hennawi & Spergel (2005). 
They multiplied (17) with a Gaussian of certain scale ra- 
dius in order to suppress the effects of the cosmic shear 
that become dominant for larger radii. Even though the 
mathematical descriptions are different, the latter two fil- 
ters are in effect very similar to (15), and we could not 
find one of them superiour over the other based on our 
rather inhomogeneous data set. Hence, we do not consider 
them for the rest of the work. The validity of our approach 
has recently been confirmed by Maturi et al. (2006), who 
also use a filter following the tangential shear profile, in- 
dividually adapted for each field to minimise the effect of 
cosmic shear. Based on the same data as we use in the 
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Fig. 8. Expected optimal S/N ratios for NFW dark matter 
haloes for four different cluster masses (M200) and two differ- 
ent image depths. The mathematical derivation of the S/N for 
a particular cluster at a given redshift is given in Appendix 
A. Note that the filter scale is not constant along each of the 
curves since the shear fields get smaller in angular size with 
increasing lens redshift. Note also that we used a lower inte- 
gration limit of z d = 0.2 for the lensed background galaxies 
(see Sect. 3.5 for both aspects). 

present paper, they find that our filter defined in equation 
15 yields very similar results as compared to their opti- 
mised filter, which means that the lensing effects of the 
large scale structure in our rather shallow survey arc not 
very dominant. 

Differences in the efficiency of such "tangential" fil- 
ters are thus expected to arise for very deep surveys only, 
and/or in case of high redshift clusters (z = 0.6 and more, 
for which our survey is not sensitive). In all other cases 
they are hardly distinguishable from each other since the 
noise in the images and the deviations from the assumed 
radial symmetry of the shear field and the NFW profile 
are dominant. Thus we consider the Qtanh filter to be 
optimally suited for our survey. For a comparison with 
other filters that do not follow the tangential profile, see 
the example shown in Fig. B.3. 

3.5. Sensitivity 

Figure 8 gives an idea of the sensitivity of our selection 
method for the Qtanh filter (see Appendix A for mathe- 
matical details). This plot was calculated taking into ac- 
count various characteristics of our survey and analysis. 
First, we introduced a maximum aperture radius of 20' to 
reflect the finite field of view of our fields. This affects 
(lowers) the S/N of very low redshift (z < 0.09) clusters, 
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whose shear fields can be very extended, and leads to a 
more distinct maximum of the curves. Second, the aper- 
ture size is not a constant along these curves but limited 
by the angular distance where the strength of the tangen- 
tial shear drops below 1%, merging into the cosmic shear. 
Lastly, we used a lower redshift cut-off of Zd = 0.2 in order 
to roughly reflect our selection criteria for the background 
galaxies (see Fig. 6). The effect of the latter cut is minor, 
it increases the S/N on the order of 5% for the low redshift 
range we probe, as compared to no cut. 

As a result, our S-statistics is insensitive for structures 
of mass equal or less than M 2 oo = 5 x 10 13 M Q for all 
redshifts. In data of average depth (n = 12 arcmin~ 2 ) 
we can detect mass concentractions of 1,2,4 x 10 14 M Q 
out to z = 0.10,0.22 and 0.32, respectively. The same 
objects would still be seen at z = 0.22,0.34 and 0.46 in 
the deeper exposures with twice the number density of 
usable background galaxies. 

3.6. Introducing the P -statistics 

In addition to the 5-statistics defined above, we intro- 
duce a new measure which we call peak position probabil- 
ity statistics, hereafter simply referred to as P-statistics. 
It tells us if at one particular position on the sky the S- 
statistics makes significantly more detections for various 
filter scales than one would expect in the absence of a lens- 
ing signal. As we find in Sect. 4, the S- and P-statistics 
complement each other rather well. 

We calculate the P-statistics as follows: 

— For each given survey field we look up the positions of 
all peaks detected in the full (6,x c ) parameter space. 

— We then select all positive peaks with a signal-to-noise 
larger than 2.5, and map their positions on the sky. The 
individual positions are weighted with the correspond- 
ing peak S/N. This discrete map is then smoothed 
with a Gaussian kernel of width 40" , yielding a con- 
tinuous map. The smoothing is introduced to smooth 
out the variations in the positions determined for the 
very same peak using different filter scales. In other 
words, it links two detections made nearby on the sky 
in different filters to the same mass concentration. 

— We randomise the orientations of the galaxies in the 
original catalogue 10 times and repeat the above step. 
From that we derive a constant noise level for each 
survey field. The assumption of a constant noise is an 
approximation which we based on a set of 200 ran- 
domisations of an arbitrary survey field. We find that 
the noise is constant to within 15% or less. It can be 
well estimated with a set of just 10 randomisations. 

— The "true" map is then divided by the rms obtained 
from the noise map, yielding an estimate of how re- 
liable the peak at the given position is. We then call 
this rescaled map the P-map in analogy to the S-map. 
However, note that the P-statistics has a very non- 
Gaussian probability distribution function (see also the 
left panel of Fig. 9). 



The main idea behind this approach is that a real peak 
has an extended shear field, i.e. it will be picked up by a 
larger number of different filter scales. In other words, as 
the aperture size changes, different samples of galaxies are 
used and all of them will yield a signal above the detection 
threshold (provided a sufficient lensing strength). On the 
contrary, a spurious peak mimicked by the noise of the 
intrinsic galaxy ellipticity is not expected to show such 
a behaviour, thus the P-statistics will prefer a true peak 
over a false peak. It has an advantage over the 5-statistics 
since it looks at a broader range of filter scales instead of 
one single scale. Thus it is capable of giving significance to 
a peak that otherwise goes unnoticed by the S'-statistics. 
On the other hand, objects with weak shear fields will not 
be recognised by the P-statistics since they appear only 
for one or a few nearby filter scales. 

We would like to emphasise that we introduced the P- 
statistics for this work on an experimental basis only. Its 
performance has not yet been evaluated based on simula- 
tions, but it yields very similar results as the 5-statistics 
(see Sect. 4), so that we included it in this presentation. 

There is some arbitrariness in the way we implemented 
the P-statistics. For example, the lower threshold of 2.5cr 
for the peaks considered can be decreased or increased. 
The former would make it smoother since more peaks are 
included, but does not yield any further advantage since it 
picks up too much noise. Increasing the threshold beyond 
3.5 reduces the number of peaks entering the statistics sig- 
nificantly. This makes the determination of the noise level 
unstable, and one starts losing less significant peaks. There 
is also room for optimisation concerning the selection of 
input data, as for this work we simply included peaks from 
the full parameter space. Concentrating on a smaller set of 
filter scales could yield a more discriminative power, but 
bears the risk of losing objects. In addition, the smooth- 
ing length has been chosen to obtain the best compromise 
between smoothing out position variations in the lensing 
detections while maintaining the spatial resolving power 
of the P-statistics. The chosen kernel width of 40" appears 
optimal for our survey, but may well be different for other 
data sets. 

In order to distinguish between individual S/N mea- 
surements made with the P- and 5-statistics, we will use 
the terms v v and v s henceforth. 

3.7. Validation of the S- and P-statistics 

As a consistency check for the concept of the 5- and the 
P-statistics, one can compare the peak probability distri- 
bution function (PDF) of the observed data against ran- 
domised data sets. The presence of cluster lensing should 
distort the PDF in the sense that more peaks are detected 
for higher S/N values (see Miyazaki et al. 2002, for an ex- 
ample) . 

To this end we created 10 copies of our entire sur- 
vey catalogue with randomised galaxy orientations, de- 
stroying any lensing signal, but keeping the galaxy po- 
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Fig. 9. Left panel: The normalised PDF for the peaks based on the entire survey area, averaged over all scales and x c . Whereas 
the peak PDF for the 5-statistics can be well approximated by a Gaussian, the PDF for the P-statistics is highly non-Gaussian 
with a broad tail. The v v should therefore not be directly interpreted as S/N. The middle (S'-statistics) and right panel (P- 
statistics) show the difference between the observed peak PDF and the PDF obtained from 10 randomisations. The middle plot 
is magnified by xlOO for \u a \ > 4 for better visualisation. Both these plots also contain the minimum (and mostly negative) 
peaks, i.e. underdense regions or voids, which accounts for the symmetric appearance. A significant excess of peaks and voids 
exists in the observed data (as compared to the randomised data) for \v a \ > 2 and \v p \ > 3. The error bars are taken from the 
randomisations. 



sitions and thus all other data characteristics fixed. We 
then calculated the PDFs for all local maxima (overdense 
regions) and minima (underdense regions) of the observa- 
tions and the randomisations, accumulating the detections 
from the entire parameter space probed. The middle and 
right panel of Fig. 9 show the difference between the PDFs 
of the observed and the randomised data sets. Both PDFs, 
for the S- and the P-statistics, are significantly skewed, 
showing an excess of peaks and voids above thresholds of 
about v s > 2 and v v > 3. 

4. Shear-selected mass concentrations 

4.1. Selecting the cluster candidates 

The way we established our sample of possible mass con- 
centrations ( "peaks" ) is as follows. The ^-statistics is eval- 
uated 

— on a grid with a 10" spacing 

— for 19 different filter scales 4 . 

— for 7 different scale parameters x c 5 

— From these S'-maps (133 different representations per 
survey field) we include all peaks higher than 4(7, which 
corresponds to the lower limit of the (S/N) m - m = 4 ... 5 
range that is considered useful by various authors (see 
Reblinsky & Bartclmann 1999; Hamana et al. 2004, for 
example) . 

The P-statistics is evaluated on the same grid space 
and obtained as described in Sect. 3.6. We will use as well 

4 1.6, 2.0, 2.4, 2.8, 3.2, 3.6, 4.0, 4.8, 5.6, 6.3, 7.1, 7.9, 8.7, 9.9, 
11.9, 13.9, 15.9, 17.9, 19.8 arcminutes. The odd numbers arise 
from the original definition which was made in units of pixels. 

5 0.025, 0.05, 0.1, 0.2, 0.5, 1.0, 2.0 



a detection threshold of v. p = 4. Oct, which has shown to 
provide us with a similar number of detections as the S- 
statistics, so we can compare the two samples. 

The detections made with both statistics are sum- 
marised in Tables B.2 to B.4. Those peaks seen with the 
S'-statistics have the best matching filter scale reported, 
i.e. the one yielding the highest S/N. Column 1 contains 
numeric labels for the detections, followed by a string in- 
dicating with which statistics it was made. The third and 
fourth column contain the detection significances f p and 
v s (if applicaple). Columns five, six and seven carry a clas- 
sification parameter (see below), a richness estimate of a 
possible optical counterpart if present, plus the distance 
of the peak from the latter. This is followed by the filter 
scale and the x c parameter (if applicable), and then the 
name of the survey field in which the detection was made. 
Finally, we report the redshift of an optical counterpart if 
known. 



4.2. Classification of the peaks detected 

In order to characterise the line of sight for a peak in 
terms of visible matter, we introduce a rough classifica- 
tion scheme based on the P-band images. We visually in- 
spected a radius of 2'.0 around the peak position for ap- 
parent overdensities of brighter galaxies, as compared to 
the surrounding field. Even though this is a rather crude 
approach, it is good enough to tell if a peak is likely asso- 
ciated with some luminous matter, or not. The radius of 
2'.0 has been chosen since we observed from known galaxy 
clusters that the lensing detection can scatter up to 1'.5 
with respect to the optical center of the cluster (see Table 
1). This is either due to substructure in the cluster, or 
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due to noise (the detection is made from a hnitc number 
of galaxies only) . 

The classes are defined as follows and are based on 
galaxies taken from the range R ~ 17 — 22 (thus they are 
not member of the catalogue of background galaxies, see 
also the left panel of Fig. 6): 

— class 1: concentration of more than 50 galaxies 

— class 2: 35-50 galaxies 

— class 3: 25-35 galaxies 

— class 4: 15-25 galaxies 

— class 5: 5-15 galaxies 

— class 6: no counterpart visible 

Judging from our number counts (Hildebrandt et al. 
2006) for the ESO Deep Public Survey data, we expect a 
number density of about 11000 galaxies with R < 22 per 
square degree. Assuming a random distribution of these 
foreground galaxies and neglecting clustering effects, we 
obtain a scattering of a = 7.2 for the total number of 
galaxies within the 2 arcmin radius. Thus, class 5 objects 
represent not more than a 2.1ct overdensity as compared 
to the randomised distribution. 

We consider classes 1 to 4 to be reliable optical counter- 
parts, and refer to them as bright peaks henceforth. Classes 
5 are rather dubious, and go as dark peaks together with 
those lines of sight classified as 6. See Fig. B.5 for an il- 
lustration of bright peaks of classes 1 — 4. 

The boundaries between the classes are permeable. For 
example, if we find an overdensity of 12 galaxies, and 4 or 5 
of them stand out from the rest by their brightness and are 
of elliptical type, the class 5 object would become a class 
4. Similarly, if we find 20 galaxies of similar brightness, 
but they show a significantly higher concentration than 
the rest of the sample, it becomes class 3. On the other 
hand, if the distance between the mass peak and the center 
of the optical peak exceeds 100 arcseconds, we decrease 
the class by one step. The same holds if the galaxies seen 
appear to be at redshift higher than ~ 0.3 or more, then 
we lower the rank by 1 since our selection method becomes 
less sensitive with increasing redshift. About 20% of our 
sample were up- or downgraded in this way. 

4.3. Spectroscopically "confirmed" candidates 

For 22 out of the 72 bright peaks, we found spectra in 
the NASA Extragalactic Database 6 . Most of them come 
from the SDSS 7 or the Las Campanas Distant Cluster 
Survey (Gonzalez et al. 2001). For 7 cases we have only 
two spectra, thus they just indicate a cluster or group 
nature, but we can not take it as hard evidence. Those 
redshifts are marked with an asterisk in Tables B.2 to B.4. 
For 15 other peaks the cluster nature was already known 
or has meanwhile been secured, either by spectroscopy, 
photometric redshifts, or by other photometric means (see 
e.g. Gladders & Yee 2000, for the red sequence method). 

6 http://ned.ipac.caltech.edu 

7 http://www.sdss.org 



Table 1. Average angular offsets between the peak and the 
optical counterpart 



class 


number 
of peaks 


average 
offset [ " ] 


1 


3 


53 ±38 


2 


14 


59 ±24 


3 


21 


51 ± 16 


4 


34 


55 ±30 


5 


43 


62 ± 32 


6 


43 


N/A 



Three of them (#043, 053 and 157) turn out to be clusters 
or groups of galaxies at different redshifts projected on 
top of each other, with the previous two being triple. For 
simplicity, we refer to all these objects in the following as 
"confirmed" peaks, even though more spectroscopic data 
has to be obtained for most of them for sufficient evidence. 

Whenever spectra were available, they confirmed our 
assumption of spatial concentrations in all cases. In order 
to secure the 50 most promising candidates, we recently 
started a large spectroscopic survey aiming at between 20 
and 50 galaxies per target. This will not only pin down the 
redshifts of the possible clusters, but also allow us to iden- 
tify further projection effects and in some cases possible 
physical connections with nearby peaks (e.g. #029, 056, 
074, 084, 092, 128, 136, 141, 158). We will report these 
results in future papers. 

The lcr redshift range of the peaks confirmed so far is 
z = 0.09. ..0.31. We therefore predominantly probe the 
lower redshift range of clusters, which is consistent with 
the theoretically expected sensitivity of our survey (see 
Kruse & Schneider 1999, and Fig. 8). 

4.4. Positional offsets 

Table 1 shows that the peaks coincide with the positions 
of the optical counterparts to within 0'.9±0'.5, indepen- 
dent of the peak classification. On the one hand, these 
offsets are due to noise, since the shear field is obtained 
by a finite number of lensed galaxies with intrinsic elliptic- 
ities. In addition, in general the shear fields deviate from 
the radial symmetry assumed by the M ap filter. On the 
other hand, these offsets can be physical in the sense that 
light does not trace mass very well for young or still non- 
virialised clusters. The two largest clusters in our sample, 
Abell 901 (#039) and Abell 1364 (#082), are good exam- 
ples. For them, the positions of the weak lensing detections 
are shifted away from the cD galaxies in the direction of 
sub-clumps. See also Fig. B.5 for an illustration. 

4.5. Comparing the S- and P-statistics 

We make 90 and 95 detections with the S- and the P- 
statistics, respectively, having 26 peaks in common. As 
can be seen from Table B.6, the number of bright and 
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Fig. 10. Shown is the spatial distribution of the 158 detected 
peaks for the WFI field of view. Open symbols indicate the 
bright peaks with classes 1-4, the filled ones peaks with classes 
5 and 6. The symbol size represents the detection significance. 
The pattern is indistuinguishable from a random distribution, 
and we also do not see differences for peaks obtained with 
either the S- or the P-statistics (not shown). 



dark peaks is equally balanced between the two meth- 
ods, both yielding slightly more dark than bright peaks. 
A few fields exist in which either only the S-statistics (e.g. 
CL1216-1201) or the P-statistics (e.g. FIELD17.P3) make 
detections, but the individual number of detections made 
per field are in general small. The statistics do not show 
a preference with respect to particular fields. The same 
holds for the spatial distribution of peaks as a function of 
position in the detector mosaic (see also Fig. 10), or the 
occurrence of dark peaks as a function of exposure time 
(see Fig. 12 for a merged plot of the two statistics). In 
general, the overlap between the two statistics increases if 
the detection thresholds of S/N — 4 are further lowered, 
i.e. peaks seen in only one of the statistics then appear 
also in the other. Yet then the overall number of detec- 
tions increases to many hundreds, and we did not evaluate 
the overlapping fraction in this case. 

The only noteworth difference between the S- and the 
P-statistics is that the latter detects about 30% peaks 
of classes 2 and 3, whereas the S'-statistics returns 25% 
more peaks of class 4. The occurence of dark peaks is 
indistuinguishable between the two methods (see Fig. 11). 

4.6. Detector and survey field biases, dark peaks 

To check for systematics concerning the WFI@2.2 detector 
array, we plotted the positions of all peaks with respect to 
the array geometry (Fig. 10). It appears that the right half 



Fig. 11. Shows the fraction of detections made in either the 
S- or the P-statistics, or in both. 



is a bit more crowded than the left part of the detector 
array. After running a dozen random distributions with 
the same number of objects we find that this is insignifi- 
cant. Thus the distribution is random- like for both bright 
and dark peaks, and does not prefer or avoid particular 
regions. 

Upon counting the bright and dark peaks in the five 
main data sources of our survey (see Sect. 2.1), we find dif- 
ferences in the ratio between bright and dark peaks (Table 
B.5). Namely, the ASTRO VIRTEL and EIS data, and our 
own observations, show an excess of 20 — 50% in terms 
of dark peaks as compared to the bright peaks, and are 
about comparable to each other. The EDisCS survey has 
a factor 2.1 more dark peaks, but is also that part of our 
survey with the most shallow exposure times. Contrary, 
the COMBO- 17 data has twice as many bright as darks 
peaks, but this comes not as a surprise since the Sll and 
A901 fields are centered on known galaxy clusters with 
significant sub-structure. If we subtract the known clus- 
ters and all detections likely associated with them, we still 
have an "excess" of 40% for the bright peaks in COMBO- 
17. Again, this is not unplausible since the COMBO-17 
fields form by far the deepest part of our survey, which 
let us detect more mass concentrations. But this holds for 
both bright and dark peaks, as the number of detections 
per square degree shows (Table B.5). 

In order to check if the dark peaks might arise due to 
imperfect PSF correction, we compare their occurrence 
with the remaining PSF residuals in our lensing cata- 
logues (Figs. 4 and 5). Therein we do not find evidence 
that the imperfect PSF correction gives rise to dark peaks. 
However, Fig. 12 indicates that small exposure times (less 
than 10 — 12 ksec) and/or a low number density (less than 
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Fig. 12. Shown is the number of bright peaks (solid line) and 
dark peaks (dashed line) as a function of exposure time (left) 
and galaxy number density (right). Shallow exposures with low 
number density have more dark peaks than deeper exposures. 



Fig. 13. Left: Histogram of the exposure times. The peak at 
57 ksec represents the Chandra Deep Field South (CDF-S). 
Right: Number density of galaxies in the 58 fields after all 
filtering steps, leaving a total number of about 710000 usable 
galaxies. The distribution reflects the distribution of exposure 
times shown in the left panel. 



n ~ 13—15 arcmin~ 2 ) of galaxies foster the occurrence 
of dark peaks. Yet this has to be seen with some cau- 
tion, in particular because we have only a small number 
of deep fields (mainly COMBO- 17) as compared to the 
shallow ones, and the deep fields are prtially concentrated 
on known structures. If we take the known structures into 
account and remove them from the statistics, we are still 
left with a smaller fraction of dark peaks in the deep expo- 
sures, but the question remains in how far the particular 
pointings of those fields introduce a bias. To answer this 
question empirically, we would need about 10 empty fields 
of 20 ksec exposure time each. 

Due to the M ap filter we use, and to the large inho- 
mogeneity of our survey, we can not directly compare the 
occurrence of dark peaks in our data to existing numeri- 
cal simulations. Also, these simulations usually make sig- 
nificantly more optimistic assumptions in terms of usable 
number density of galaxies and field of view than we could 
realise with GaBoDS (see Reblinsky & Bartclmann 1999; 
Jain & van Waerbeke 2000; Hennawi & Spergel 2005, for 
example). In particular Hamana et al. (2004) have shown 
that in their simulations (n = 30) they expect to detect 43 
real peaks (efficiency of <~ 60%), scaled to the same area 
as GaBoDS and drawn with S/N > 4 from mass recon- 
struction maps. A similar number of false peaks appear 
as well, being either pure noise peaks, or peaks with an 
expected S/N < 4 being pushed over this detection limit. 
The latter would be labelled as bright peaks in our case. 
Our absolute numbers are different (72 bright and 86 dark 
peaks) since we use a very different selection method. Yet, 
if we interpret our dark peaks as noise peaks, the ratio be- 
tween our bright and dark peaks is comparable to the ratio 
between their true and false peaks. 

This interpretation, i.e. dark peaks are mostly noise 
peaks, is strengthened by the fact that with increasing 
peak S/N the fraction of dark peaks is decreasing (see 
Table B.6), for both the S- and the P-statistics. However, 
our observational data base (sky coverage) is too small to 
tell if this trend, i.e. the dark peaks dying out, continues 
for higher values of the S/N. 



5. Summary 

In the present paper we have introduced a new type of 
filter function for the aperture mass statistics. This fil- 
ter function follows the tangential shear profile created by 
a radial symmetric NFW density profile. We have shown 
that it is optimally suited for an application to our 19 
square degree weak lensing survey conducted with the 
WFI@2.2m MPG/ESO telescope. This is a survey with 
very inhomogeneous depth, in which we expect to be 
able to detect mass concentrations in the redshift range 
z = 0.1 . . .0.5. We also introduced the new P-statistics, 
currently on a test-basis only. It turned out to deliver very 
similar results in the number of bright and dark peaks de- 
tected, and for the time being looks like a good comple- 
ment to the S'-statistics. Its performance has to be inves- 
tigated more deeply though, and there is space for optimi- 
sation. The global PDFs for both the S- and P-statistics 
show clear excess peaks for higher values of S/N as com- 
pared to the randomisations. Thus the presence of lensing 
mass concentrations in our survey data is confirmed. 

We introduced a classification scheme in order to asso- 
ciate the hypothetical mass peaks detected with possible 
luminous matter. From the 158 detections we made with 
the combined S- and P-statistics, 72 (46%) appear to have 
an optical counterpart. For 22 of those we found spectra 
in the literature, confirming the above mentioned redshift 
range, and that indeed a mass concentration exists along 
those particular lines of sight. For a smaller number of 
those we have spectroscopic evidence that the peaks de- 
tected are due to projection effects. We expect that in our 
currently conducted spectroscopic follow-up survey more 
such projection cases will be uncovered, together with a 
confirmation of a very significant fraction of the remaining 
bright peaks. In a future paper we will also compare this 
shear-selected sample with an optically selected sample 
using matched filter techniques. 

We gained some insight into the subject of dark peaks. 
They appear preferentially in shallow data with a small 
number density of galaxies, indicating that a large frac- 
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tion of those could be due to noise (i.e. instrinsic galaxy 
cllipticies), or they are statistical flukes. Nevertheless, in 
our deep fields we also observe a significant fraction of dark 
peaks, but our statistics can not be unambiguously inter- 
preted since those fields are biased towards clusters with 
significant sub-structure, and we have only a very small 
number of them. Nevertheless, real physical objects such 
as very underluminous clusters are far from being ruled 
out at this point. At least on the mass scale of galaxies 
the last year has seen astonishing examples of objects hav- 
ing several 10 8 Mq of neutral hydrogen, yet apparently 
no star formation has taken place in them (Minchin et al. 
2005; Auld et al. 2006). Whether similar objects can still 
exist on the cluster mass scale is unclear. 

Finally, we would like to repeat that the Garching- 
Bonn Deep Survey has been made with a 2m telescope. 
Most numerical simulations done so far are much more 
optimistic in terms of the number density reached (n ~ 
30 arcmin~ 2 ) and correspond to surveys that are cur- 
rently conducted (or will be in the near future) with 4m- 
and 8m-class telescopes, such as SUPRIME33 (Miyazaki 
et al. 2005) or the CFHTLS 8 . One noteworthy exception 
will be the KIDS 9 (Kilo Degree Survey) obtained with 
OmcgaCAM@VST, covering 1500 square degrees starting 
in 2007. 
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Appendix A: Expected S/N for a NFW halo 

We derive the S/N expected for a radially symmetric 
NFW dark matter halo, using a flat cosmology with 
n = 0.3, n x = 0.7, h = 0.7 and er 8 = 0.85. 

According to S96, the S/N for M ap for a cluster at the 
origin of the coordinate system evaluates as 



s{0) = /47m / P d^T? 7 tWQW 

V ^ /]>W 



(A.l) 



Here we assumed radial symmetry (yielding a factor of 
2-7r) , write 9 for the aperture radius and •& for the distance 
from the centre of the aperture, n is the number density 
of galaxies, and a £ is the dispersion of the modulus of the 
galaxy ellipticities. 

The tangential shear for a radially symmetric NFW 
profile was given by Wright & Brainerd (2000) as 
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Here, cl is the speed of light and G is Newton's con- 
stant. The NFW concentration parameter c is a function 
of cosmology, the mass of the cluster and its redshift. We 
calculate c using the cens routine kindly provided by J. 



Navarro 



10 



and 5 r are the NFW scale radius and charac- 



teristic over-density of the cluster. S cr is the critical sur- 
face mass density, where we write Dd, D s and D^s for the 
angular diameter distances between the observer and the 
lens, between the observer and the source, and between 
the lens and the source, respectively. For a flat cosmology, 
they are defined as 



D(z 1 ,z 2 ) = — — / da [a fl + a 4 fi A ] 

-no 1 + Z 2 J a2 



-1/2 



(A.7) 

with z\ < z 2 and a — 1/(1 + z). To take into account the 
redshift distribution of the lensed galaxies, we assume that 
those follow the normalised distribution given in Brainerd 
et al. (1996), with parameterisation 



^ ^ 2zq (zq 



exp 



3/2' 



10 http://www.astro.uvic.ca/~jfn/cens/ 



(A.8) 



The ratio D ds /D B is averaged over this redshift distribu- 
tion, starting with the lens redshift lower integra- 
tion limit. The latter was chosen because galaxies with 
z < Zd are unlensed and largely removed from our cata- 
logues by appropriate detection thresholds and cut-offs 11 . 

The functional expression for g(x) is identical to the 
one already given in equation (14), and contains the shape 
of the shear profile. Finally, fixing the remaining numerical 
parameters provides us with all information to calculate 
the S/N. From our data we have a e = 0.48, and we assume 
two different image depths which we base on empirical 
findings. One is shallow with n = 12 arcmin~ 2 and z = 



0.4, and the deeper one given by n 
zo = 0.5. 

The S/N then evaluates as 



24 arcmin and 



V a * Scr ^// o 9 di?0Q 2 (0) 



(A.9) 



11 The exact effect of this filtering on the shape of the redshift 
distribution is not yet investigated. It will mostly affect the 
S/N prediction of high-redshift clusters, and only very little 
the low-z regime probed with our data. 
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Appendix B: Further tables and figures 



Table B.l. The 58 fields used for this work 



Field 


a(2000.0) 


5(2000.0) 


T oxp 


Source 


A1347.P1 


175.257 


-25.514 


13500 


own observation 


A1347.P2 


175.792 


-25.509 


7500 


own observation 


A1347.P3 


175.239 


-25.009 


7000 


own observation 


A1347.P4 


175.794 


-24.998 


8000 


own observation 


A901 


149.077 


-10.027 


18100 


COMBO-17 


AMI 


58.811 


-49.667 


7500 


ASTROVIRTEL 


B8ml 


340.348 


-10.089 


4500 


ASTROVIR.TEL 


B8m2 


340.348 


-10.589 


5400 


ASTROVIRTEL 


P»8m3 


340.346 


— 11.088 


5400 
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B8p0 


340.348 


-9.590 


7200 
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B8pl 


340.346 


— 9.089 


4500 
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B8p2 


340.345 


-8.589 


5400 
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B8p3 


340.345 


-8.089 


5400 
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214.360 


-12.253 


4800 


ASTROVIRTEL 


C04ml 


214.727 


-12.753 


4000 


ASTROVIRTEL 


C04m2 


214.478 


— 13.253 


4000 


ASTROVIRTEL 


C04m3 


215.318 


— 13.753 


4000 


ASTROVIRTEL 


C04m4 


215.111 


-14.253 


4000 


ASTROVIRTEL 


C04pl 


214.726 


-11.753 


4000 


ASTROVIRTEL 


C04p2 


214.727 


-11.253 


4000 


ASTROVIRTEL 


C04p3 


215.098 


— 10.753 


4000 


ASTROVIRTEL 


CAPO-DF 


186.037 


— 13.107 


13000 


ASTROVIRTEL 


CDF-S 


53.133 


— 27.822 


57200 


GOODS EIS 
COMBO-17 


CL1037— 1243 


159.444 


— 12.754 


3600 


EDisCS 


CL1040— 1155 


160.139 


— 11.963 


3600 


EDisCS 


CL1054— 1146 


163.581 


— 11.813 


3600 


EDisCS 


CL1054-1245 
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-12.797 


3600 
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CL1059 — 1253 


164.755 


-12.920 


3000 
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-11.525 


3600 


EDisCS 
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174.508 


— 11.599 
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180.645 


-12.441 


3600 


EDisCS 
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184.170 
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— 11.630 
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208.306 
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215.066 
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52.506 


-27.817 
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DEEP2f 


53.669 


-27.324 
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171.245 
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-34.694 


10000 


own observation 


F17.P3 
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-1.734 
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11.498 


-29.610 
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76.333 


-28.818 


11400 
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Fig. B.l. Indicated are the objects that are left over in the lensing catalogue after all typical filtering steps. Only fainter 
background galaxies are kept. Brighter sources, spurious detections, stars and highly elliptical objects such as asteroid tracks 
are largely absent from the catalogue. The field of view is about 3' . 




Fig. B.2. Typical PSF anisotropy. Upper row: £1,2 scatter plot for stars before and after PSF correction. Lower left: PSF 
anisotropy pattern before correction. This plot is a more intuitive representation of the left panel above. Lower right: PSF 
residuals after a polynomial fit was used to correct for the anisotropy. No coherent shear signal is left. 
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Fig. B.3. Tangential shear and S'-profile for the two largest clusters in the survey data. In the upper row the tangential shear 
is shown, with the TANH filter that yielded the highest S/N overlaid as a solid line. For better comparison the amplitude of 
TANH was scaled so that it best fits the tangential shear. All data points are mutually independent. The lower row shows the 
S-profile of the two clusters for different filters Q. The NFW filter is plotted for 10 different values of x c € [0.01, 1.5], to show 
the scatter delivered by x c . The filters defined by Padmanabhan et al. (2003) and Hennawi & Spergel (2005) (not shown) deliver 
very similar results. For comparison we plot various other types of filter functions introduced in the literature. POLY is the 
polynomial filter defined by Schneider et al. (1998), mainly as a new measure for cosmic shear rather than cluster detection. 
The filters S96 were defined by Schneider (1996). EXP is based on the difference of two Gaussians of different width (Schirmer 
2004). Clearly, all of them yield smaller S/N values than those following the tangential shear profile. If one of those filters yields 
a higher S/N for a cluster than TANH in our data, then this is marked accordingly in Tables B.2 and B.3. 
Note: Although the tangential shear is smaller for A901 than for Sll, the S/N is higher due to the larger number density of 
galaxies with measured shapes (n = 15 arcmin~ 2 for Sll, and n — 24 arcmin~ 2 for A901). 
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Fig. B.4. Upper right: S-map for randomised galaxy orientations and positions. Upper left: The same S-map, evaluated after 
10 randomly positioned holes with a radius of 90% of the filter scale were cut into the data field (lower left). Artificial features 
show up in the S-map at the positions of the holes, the latter being exaggerated in number and size for better visualisation. 
Lower right: True galaxy distribution of the field from which the galaxies were drawn for this test. The largest hole is caused 
by an 8th magnitude star. We conclude that holes in the data fields are in general not a cause of concern for our analysis. 
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Table B.2. Shear-selected mass concentrations (part 1). The first column contains a label for the peak, and the second 
one indicates if the detection was made with either the S- or the P-statistics (or both). The next two columns contain the 
corresponding significances. The classification shows if along the line of sight an overdensity of galaxies is found, with class 
1 meaning a very obvious overdensity, and class 6 no overdensity. The richness indicates how many galaxies were found as 
compared to the average density in this field, and we give the distance between the peak and the optical counterpart. Finally, 
we have the name of the survey field in which the detection was made, and possibly a redshift for the counterpart. An asterisk 
behind the redshifts indicates that it is based on the measurement of less than three galaxies. 
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3 


15 


50 


6'.3 


2.0 


A901 




#036 


p 




7.08 


2 


45 


60 






A901 


0.161 


#037 


p 




5.01 


4 


15 


60 






A901 




#038 


p 




4.60 


4 


25 


50 






A901 




#039 


SP 


6.30 


10.01 


1 


50 


70 


19'.8 


0.1 


A901 


0.161 


#040 


p 




6.08 


3 


35 


50 






A901 


0.165 


#041 


p 




4.52 


4 


20 


80 






Pal3 




#042 


s 


4.13 




3 


25 


40 


2'.8 


2.0 


Pal3 


0.189* 


#043 


p 




5.14 


4 


20 


40 






Pal3 


0.094, 0.123, 0.215 


#044 


p 




4.26 


3 


30 


40 






CL1037-1243 




#045 


p 




5.00 


6 










CL1040-1155 




#046 


s 


4.19 




2 


40 


70 


2'.4 


0.5 


CL1040-1155 




#047 


p 


3.53 


4.06 


5 


5 


120 






CL1040-1155 




#048 


s 


4.07 




6 






9'.9 


1.0 


CL1040-1155 




#049 


s 


4.18 




5 


5 


80 


5'.6 


0.1 


CL1054-1245 


0.122* 


#050 


s 


4.07 




6 






3'.6 


0.05 


CL1054-1146 
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Table B.3. Shear-selected mass concentrations (part 2). 











ClctSS 


rich 
nc 


LLlbLr 


sciilc 


x c 


FiplH 




#051 


S 


4.21 




5 


15 


20 


5'.6 


0.1 


CL1059-1253 




#052 


p 




4.60 


5 


20 


30 






CL1059-1253 




#053 


s 


4.02 




4 


20 


120 


6'.3 


2.0 


CL1059-1253 


0.41, 0.52, 0.60 


#054 


s 


4.01 




4 


15 


30 


2'.4 


2.0 


CL1059-1253 




#055 


SP 


4.53 


6.31 


5 


5 


20 


3'.2 


2.0 


DEEP3D 




#056 


p 




4.08 


4 


15 


10 






DEEP3D 




#057 


SP 


4.24 


7.48 


4 


25 


80 


5'.6 


0.1 


DEEP3D 




#058 


SP 


4.38 


5.04 


4 


20 


70 


5'.6 


0.5 


CL1119-1129 




#059 


s 


4.25 




6 






1'.6 


0.1 


DEEP3D 




#060 


p 




5.82 


6 










DEEP3D 




#061 


s 


4.08 




4 


20 


40 


2^0 


2.0 


DEEP3A 




#062 


s 


4.21 




4 


20 


70 


3'.2 


0.2 


DEEP3A 




#063 


p 




4.67 


5 


20 


120 






CL1138-1133 




#064 


p 




4.31 


4 


15 


60 






CL1138-1133 




#065 


p 




4.97 


4 


30 


50 






CL1138-1133 




#066 


SP 


4.02 


5.48 


3 


25 


50 


7'.9 


0.05 


A1347.P1 




#067 


p 




4.10 


5 


35 


50 






A1347.P3 




#068 


s 


4.00 




5 


5 


50 


1'.6 


0.1 


A1347_P1 




#069 


SP 


4.21 


4.95 


5 


10 


60 


6'.3 


0.025 


A1347.P3 




#070 


SP 


4.32 


4.94 


6 






7'.1 


0.05 


A1347.P3 




#071 


s 


4.28 




4 


20 


50 


2'.4 


2.0 


Sll 




#072 


p 




4.37 


4 


15 


90 






A1347_P1 




#073 


p 




6.17 


3 


25 


70 






A1347_P1 




#074 


s 


4.27 




4 


40 


60 


1'.6 


0.5 


Sll 




#075 


s 


4.43 




3 


35 


30 


1'.6 


0.5 


Sll 




#076 


s 


4.27 




6 






2'.8 


2.0 


Sll 




#077 


s 


4.76 




4 


20 


70 


13'. 9 


0.05 


Sll 




#078 


SP 


4.00 


4.05 


2 


40 


40 


2'.8 


2.0 


Sll 


0.119 


#079 


p 




4.02 


3 


30 


40 






A1347.P4 




#080 


s 


4.02 




4 


20 


10 


2'.8 


0.1 


A1347.P4 




#081 


p 




4.78 


6 










A1347.P2 




#082 


SP 


4.87 


5.83 


1 


60 


80 


13'. 9 


0.025 


Sll 


0.106 


#083 


p 




5.33 


5 


15 


60 






CL1202-1224 




#084 


p 




4.61 


5 


40 


120 






CL1202-1224 


0.4 


#085 


s 


4.63 




5 


10 


50 


2'.4 


0.5 


CL1216-1201 




#086 


s 


4.43 




5 


10 


90 


4'.0 


0.5 


CL1216-1201 




#087 


s 


4.01 




5 


30 


130 


5'.6 


2.0 


CL1216-1201 


0.33 


#088 


s 


4.64 




6 






3'.2 


0.2 


CL1216-1201 




#089 


p 




5.01 


5 


10 


70 






CL1301-1139 




#090 


s 


4.18 




3 


30 


70 


9'.9 


2.0 


CL1301-1139 




#091 


SP 


4.64 


4.99 


3 


30 


50 


11'.9 


0.025 


CL1353-1137 




#092 


p 




5.28 


3 


35 


50 






CL1353-1137 




#093 


SP 


4.39 


4.40 


5 


15 


90 


3'.6 


2.0 


C0400 




#094 


s 


4.14 




6 






5'.6 


2.0 


C0400 




#095 


p 




5.30 


6 










C04p2 




#096 


p 




4.32 


6 










C04pl 




#097 


p 




4.15 


6 










C04ml 




#098 


SP 


4.64 


7.60 


4 


25 


20 


5'.6 


2.0 


C04p2 




#099 


S 


4.22 




6 






2'.4 


2.0 


C04m4 




#100 


p 


5.13 


5.90 


6 










C04p3 
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Table B.4. Shear-selected mass concentrations (part 3). 



Peak 


stat 




p 


class 


rich 


dist 


scale 




Field 


z 


#101 


P 




5.59 


6 










CL1420-1236 




#102 


P 




5.59 


6 










C04m4 




#103 


P 




4.55 


6 










C04m4 




#104 


P 




4.90 


4 


15 


80 






C04m4 




#105 


S 


4.08 




6 






2'.4 


2.0 


CL1420-1236 




#106 


s 


4.07 




6 






6'.3 


2.0 


CL1420-1236 




#107 


p 




4.12 


6 










CL1420-1236 




#108 


s 


4.43 




6 






2'.4 


1.0 


C04m3 




#109 


SP 


4.03 


6.29 


5 


10 


20 


6'.3 


0.2 


CL1420-1236 




#110 


p 




4.02 


3 


35 


70 






C04p3 




#111 


s 


4.05 




5 


30 


90 


4'.0 


0.1 


CL1420-1236 




#112 


s 


4.01 




6 






4'.0 


0.05 


C04m3 




#113 


s 


4.07 




6 






2'.4 


0.5 


FIELD 17.P1 




#114 


p 




4.59 


2 


80 


60 






FIELD17_P3 




#115 


p 




4.59 


3 


25 


50 






FIELD17_P3 




#116 


s 


4.28 




6 






1'.6 


2.0 


FIELD 17.P1 




#117 


p 




4.60 


5 


10 


30 






FIELD 17_P3 




#118 


p 




4.97 


5 


10 


90 






FIELD 17.P3 




#119 


p 




4.35 


6 










FIELD17.P3 




#120 


p 




5.76 


3 


25 


40 






FIELD 17_P3 


0.132* 


#121 


s 


4.11 




5 


20 


70 


2'.0 


0.5 


FIELD4.P1 




#122 


s 


4.09 




4 


20 


40 


3'.6 


0.2 


FIELD4.P1 




#123 


s 


4.02 




3 


35 


40 


2'.0 


0.2 


FIELD4.P3 




#124 


s 


4.00 




4 


15 


10 


2'.4 


1.0 


FIELD4.P3 




#125 


p 




4.35 


6 










FIELD4.P4 




#126 


s 


4.02 




5 


5 


30 


2'.0 


2.0 


FIELD4.P3 




#127 


p 




4.46 


6 










FIELD4.P4 




#128 


SP 


4.13 


4.21 


3 


30 


60 


7'.9 


0.025 


FIELD4.P4 




#129 


s 


4.00 




5 


15 


70 


3'.2 


0.1 


FIELD4.P4 




#130 


SP 


5.10 


7.61 


4 


15 


80 


17'.8 


0.05 


B8p0 




#131 


SP 


4.29 


4.30 


5 


10 


90 


19'.8 


0.05 


B8ml 




#132 


s 


4.12 




5 


10 


70 


1A 


1.0 


B8m2 




#133 


p 




4.05 


6 










B8m3 




#134 


SP 


4.05 


6.06 


6 






7'.1 


2.0 


B8p0 




#135 


p 




4.14 


6 










B8m2 




#136 


SP 


4.45 


10.20 


3 


30 


50 


2'.8 


0.5 


B8p0 




#137 


SP 


4.12 


6.26 


6 






3'.6 


2.0 


B8p0 




#138 


p 




4.70 


5 


15 


60 






B8p0 




#139 


SP 


4.35 


6.95 


4 


15 


10 


2'.0 


0.2 


B8p0 




#140 


s 


4.13 




4 


15 


90 


15'.9 


0.05 


B8p2 


0.07* 


#141 


SP 


4.59 


7.76 


2 


40 


80 


3'.2 


0.5 


B8p0 




#142 


s 


4.71 




3 


35 


70 


4'.0 


2.0 


B8pl 


0.318 


#143 


s 


4.07 




6 






6'.3 


0.1 


B8p3 




#144 


p 




4.71 


2 


40 


80 






B8p0 




#145 


s 


4.25 




6 






1'.6 


0.1 


B8ml 




#146 


s 


4.11 




6 






3'.6 


2.0 


B8pl 




#147 


p 




6.78 


2 


90 


100 






B8p0 




#148 


p 




4.31 


4 


15 


10 






B8p0 




#149 


s 


4.27 




5 


10 





2'.0 


0.5 


B8ml 




#150 


p 




4.21 


4 


30 


80 






B8p0 




#151 


p 




4.01 


6 










B8m2 




#152 


SP 


4.14 


4.20 


5 


10 


30 


2'.0 


1.0 


B8p0 




#153 


p 




4.12 


5 


25 


100 






DEEP1C 




#154 


p 




4.33 


6 










DEEP IE 




#155 


s 


4.02 




5 


20 


120 


3'.6 


0.5 


DEEP1C 


0.135* 


#156 


p 




6.22 


6 










DEEP1A 




#157 


s 


4.35 




4 


20 


30 


7'.9 


0.05 


DEEP1A 


0.151, 0.176 


#158 


p 




6.91 


2 


40 


30 






DEEP1A 


0.175 
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Table B.5. Bright (classes 1 to 4) and dark (5-6) peaks for the various survey data sources. The columns contain the data 
source, the number of bright and dark peaks, the ratio between dark and bright peaks, the average exposure times and image 
seeing, the area covered, plus the number of bright and dark peaks per unit area. For the COMBO-17 field we give in parenthesis 
the corresponding values when the known structures are subtracted. 



Source 


n b 


n d 


n d /n h 


t [ksec] 


s["] 


A[D 2 ] 


n h /A 


n d /A 


EDisCS 


10 


21 


2.1 


3.5 


0.90 


3.8 


2.6 


5.5 


ASTRO VIRTEL 


18 


27 


1.5 


4.9 


0.87 


6.1 


3.0 


4.4 


COMBO-17 


25 (18) 


13 


0.5 (0.7) 


25.7 


0.80 


1.6 


15.6 (11.2) 


8.1 


EIS 


6 


7 


1.2 


6.8 


0.94 


2.9 


2.1 


2.4 


Own obs. 


13 


18 


1.4 


8.9 


0.88 


4.2 


3.1 


4.3 



Table B.6. Number of bright (classes 1-4) and dark (5-6) peaks and their ratios for the S- and the P-statistics. 



S/N 


n b (S) 


n d (S) 


n d /n b 


n h (P) 


MP) 


n d /n h 


4.0 - 4.25 


23 


31 


1.3 


7 


12 


1.7 


4.25 - 4.5 


9 


11 


1.2 


5 


9 


1.8 


4.5-4.75 


6 


5 


0.8 


6 


7 


1.2 


4.75 - 5.0 


2 





0.0 


3 


4 


1.3 


5.0-6.0 


2 


1 


0.5 


8 


11 


1.4 


> 6.0 


1 





0.0 


15 


7 


0.5 


total 


43 


48 




44 


50 
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Fig. B.5. Illustrates the typical appearance of bright clusters of various classes, as defined in Sect. 4.2. Note that the resolution 
is in general not high enough to distinguish smaller member galaxies from stars. The field of view is 4'.3. 



